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ABSTRACT 

We present results on the size evolution of passively evolving galaxies at 1 < z < 2 drawn from 
the Wide Field Camera 3 Early Release Science program. Our sample was constructed using an 
analog to the passive BzK selection criterion, which isolates galaxies with little or no on-going star 
formation at z > 1.5. We identify 30 galaxies in ~ 40 arcmin 2 to H < 25 mag. We supplement 
spectroscopic redshifts from the literature with photometric redshifts determined from the 15-band 
photometry from 0.22 — 8 /im. We determine effective radii from Sersic profile fits to the 7f-band 
image using an empirical PSF. We find that size evolution is a strong function of stellar mass, with 
the most massive (M* ~ 10 11 M Q ) galaxies undergoing the most rapid evolution from z ~ 2 to the 
present. Parameterizing the size evolution as (1 4- z)~ a 1 we find a tentative scaling between a and 
stellar mass ofaw —1.8 + 1.41og(A//10 9 -M©). We briefly discuss the implications of this result for 
our understanding of the dynamical evolution of the red galaxies. 

Subject headings: Keywords: galaxies: evolution — galaxies: structure — galaxies: fundamental 
parameters 



1. INTRODUCTION 



Understanding the red, passively evolving, galaxies at 
intermediate redshifts (z ~ 2) is one of the outstand- 
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ing challenges of galaxy evolution studies. Early ex- 
pectations of a high luminosity phase associated with 
the formation of spheroids on a fre e-fall time scale 
(jEggen. Lynden-Bell. fc Sandagd lT962T ) have long given 
way to a framework in which the spheroids are assem- 
bled over an extended period of time. The identifica- 
tion of fairly high space densities of massive and pas- 
sively evolving galaxies at z> 1.5 (e.g. lGlazebrook et al.1 
l2004t iCimatti et al.l I2004T ) revealed weaknesses in the 
early semi-analytic models of their formation. Large- 
scale surveys, both locally and out to z ~ 1, provide a 
fairly clear view of the evolution in the number density 
of red sequence galaxies an d the global stellar mass den- 
sity in passive systems (e.g.lBell et al.ll2006t iFaber et"aTI 
120071: IBrown et all [20031 12007ft . 

Deep spectroscopic studies of red galaxies (e.g. 
iMcCarthv et~aT1 l200l ICimatti et all 12001 l200l esti- 
mated stellar ages consistent with early formation red- 
shifts (zform ^ 4), potentially at odds with the rapid 
evolution in stellar mass density at these epochs (e.g. 
IRudnick et al.l I2003T) . It is now clear that t he red se- 
quen c e was in place at f airly early epochs (Be ll et aT] 
120041: iBremer et all 120061). but evolved strongly in the 
1 < z <3 era fe.g. iDemarco et"aT]|2010D and likely since 
z~l as well fe.g. IFaber et a l. 2007). The emerging pic- 
ture in which massive galaxies are assembled via mergers 
followed by a rapid quenching of star formation addresses 
ma ny of the salient pr operties of the red sequence galax- 
ies (|Faber et all 120071 ). 

In light of these results it was quite surprising that 
high spatial resolution studies showed that the passive 
red galaxies at intermediate redshifts are systematically 
smaller than thei r likely present-day counterparts (e.g. 
iDaddi et al.l [20051 ItTuuIIo et al.l [2001 . Observations of 
I < z < 2 passive galaxies with either Advanced Camera 
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of Surveys (ACS) or Near Infrared Camera and Multi- 
Object Spectrometer (NICMOS) on the Hubble Space 
Telescope (HST) gave sizes that are a fact or of ~ 2 smaller 
than equal mass red galaxies today (e.g.lLqnghetti et al.1 
120071 iCimatti et al.ll200l iDamjanov et all 120091 ) . while 
observations of red galaxies at z > 2 suggest even 
more dramatic evolution (e.g. Ivan Dokkum et all 120081 : 
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o: 



iSaracco. Longhetti. fc Gargiuloll2010() 

The primary interest in determining the characteristic 
sizes of passive galaxies is the insight that it provides 
into the dynamical state of t he system. As one projec- 
tion of the fun damental plane (jDjorgovski fc pavislll987l: 
Dressier et al|[l987T) . the Kormendv relation (|Kormendvl 



19771) provides a probe of the dynamics of hot stellar sys- 



tems. The Kormendy diagram for the z > 1 red sequence 
galaxies is nearly as tight as it is loca lly, but is offset 
in bo th surface brightness and size (e.g.lDamiano v et all 
2009) , even when the effects of passive fading of the stel- 
lar populations is taken into account. The compact sizes 
and high surface brightness of the red galaxies implies 
stellar densities in the centers of these galaxies that are 
two to three orders of magni tude higher than present 
day massive red galax ies (e.g. Ivan Dokkum e t al. 2008; 
IDamjanov et al.ll2009D . 

The challenge in understanding these results stems 
from the apparent conflict between the requirement for a 
factor of ~ 3 growth in size during an epoch in which the 
stellar masses, population ages, and overall morphologies 
show little or no evolution. A number of models are un- 
der active discussion, but most involve large numbers of 
late stage minor mergers that grow the galaxies in size 
without contributing large amounts of mass. While it is 
difficult to conclusively rule out various measurement bi- 
ases which might give rise to this apparent trend (e.g. un- 
derestimated effective radii or overestimated masses), 
most interpretations are astrophysical in origin: (ma- 
jor or minor) merging occurring a t different phases (e.g. 
iNaab. Khochfar. fc Burkertl I2006D . or expansion due to 
a significant mass loss, either by an ac tive galactic nu- 
cleus (|Fan et al . 2008) or stellar winds (IDamjano v et all 
2009). Recently Hopki ns et all (|2010al ) present a semi- 
analytic mo del based on hig h-resolution hydrodynamic 
simulations (|Cox et all 120061 ). which incorporates sev- 
eral astrophysical and observational mechanisms for this 
trend, and conclude that for galaxies with a stellar mass 
of M* > 10 11 Mq, the late-stage minor merging is the 
dominate mechanism, accounting for ~ 50% of the ap- 
parent size evolution. Certainly many of these high- 
reds hift systems do in fact show merger- like features 
(e.g. Ivan D okkum 2005; Ivan Dokkum fc Brammerll2010D 
or have (multiple) companions (e.g. lBell et alll2006l ). 

The merging s cenario implies that some non- negligible 
fraction (< 10%; Hop kins et alJ l2009) of c ompact galax- 
ies sh ould remain to z ~ 0. However iTrujillo et all 
(2009) find only 0.03% (= 48/152,083) galaxies at 
z ~ 0.2 have stellar densities comparable to these 
high redshift galaxies, arguing against the merger sce- 
nario. To further confound the issue, these galax- 
ies are generally young (~ 3 Gyr old), giving an ap- 
proximate formation redshift of Zf or m ~ 0.3, suggest- 
ing that they are not the relics of the early universe, 
but for med from gas-rich, recently merged dj sks. Con- 
versely, [SaraccoITonghettil^Gargiulo (2010) find that 



within la of the local R e — M* relation (She n et all 
2003). Based on this finding, they suggest that the 
compact high redshift galaxies are not the progenitors 
of large field galaxies, but of compact brightest clus- 
ter galaxies. The absence of compact early- type galax- 
ies at z ~ seems to suggest that the simple pic- 
ture of galaxy merging with some mass-loss expansion 
is not completely c orrect, in stark contrast to the sim- 
ulated results (e.g. Naab. Johansson, fc Ostrikerl 120091 : 



62% (= 



21/34) of their galaxies at 1 < z spcc < 2 are 



iHopkins et"aTll2010al )7 

This paper is organized as follows. In §[2] we describe 
the observations and ancillary data, in §[3] we define our 
sample, in §2] we detail our stellar population modeling 
and morphological analysis, in §[4] we describe our mass- 
dependent size evolution model, in §[6] we discuss several 
key results, and in § [7] we give a brief summary of this 
work with thoughts for future surveys. Throughout this 
paper we assume a ACDM cosmology wi th f2p = 0.27, 
Q A = 0.73, and H = 73 km s" 1 Mpc" 1 (|Spergel et all 
120031) and wi ll quote all magnitu des in the AB^ magni- 
tude system (|Oke fc Gunnlll983ft . 

2. OBSERVATIONS 
2.1. The Wide Field Camera 3 Early Release Science 

For this work we analyze the Early Release Science 
(ERS; PropID: 11359, PI: R. W. O'Connell) observa- 
tions conducted with the Wide-Field Camera 3 (WFC3) 
recently installed on the Hubble Space Telescope (HST). 
This field covers roughly the northern 40 arcmin 2 of the 
Great Observatori es Origins Deep Survey southern field 
(GOODS-S; Giavalisco et al. 2004), and leverages the ex- 
isting Advanced Camera for Surveys (ACS) optical data 
with near-ultraviolet (NUV) and near-infrared (NIR) 
data with equivalent space-based imaging. The main 
details regarding the data collection, r eduction, calibra- 
tion, a nd mosaicking are presented by Wind horst et all 
(|2010f ). We briefly summarize the observational aspects 
critical to this work. 

The ERS program utilizes the two complementary 
modes of the WFC3: UVIS and IR channels. The 
data were collected between September and Novem- 
ber 2009 and form a 2 x 4 and 2x5 mosaic, respec- 
tively. The final dataset consists of ~ 40 arcmin 2 of 
10-band HST imaging covering 2250 A to 1.6 /im in 
wavelength. The v2.0 GOODS high-level science prod- 
ucts had an original pixel scale of 30 milliarcseconds, 
and were 3x3 block summed to produce science and 
weight map s of equal pixel scale to the WFC3 obser- 
vations (see iWindhorst et all 120101 for a more detailed 
discussion on the rebinning process and its motivation). 
Therefore, our final HST mosaics have 90 milliarcsecond 
pixels. We refer hereafter to the 10-band WFC3 and 
ACS filter set in the F225W, F275W, F336W, F435W, 
F606W, F775W, F8501P, F098M, F125W, and F160W 
as U\UiU?,BVi' z'Y s JH to simplify the bandpass nota- 
tion. 

2.2. Source Catalogs 

We created source catalogs using SExtractor 
(jBertin fc Arnouti fl996) with the if-band mosaic as the 
detection image and each of the HST mosaics as mea- 
surement image s. We use the wei ght maps produced by 
MultiDrizzle (jKoekemoer et alj|2002f ). modified to ac- 
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count for correlated pixel noise ((Dickinson et al.l I2004D - 

For object detection, we require a minimum of 5 con- 
nected pixels each greater than 0.6a above the local 
background and apply a 3 pixel Gaussian smoothing fil- 
ter. For deblending, we adopt a contrast parameter of 
10~ 4 with 64 sub-thresholds. For photometry, we adopt 
the MAG_AUT0 measurements with a Kron factor of 2.5, 
and minimum object radius of 3.5 pix els, which reliably 
recov er total fluxes to within ~ 6% (iBertin fc Arnouts 
I1996T). We use the AB zeropoints given bv iKalirai et al. 
(2009a.b) for the UVIS and IR data, respectively. 

3. SAMPLE SELECTION 

We select our passively evolving ga laxies using a varia - 
tion on the standard pBzK selection (jDaddi et al.l l2004h 
designed to identify galaxies at z ~ 2 with little on-going 
star formation and stellar population ages of ~ 1 Gyr 
(jDaddi et al.l [2005). However, to take full advantage 
of our high-spatial resolution and significantly deeper 
ERS dataset, we use our ii-band imaging instead of the 
ground-based K-band images. Therefore we modify the 
classical pBzK selection criteria to an equivalent pBzH 
scheme: 

0' -H)-{B- z) < -0.2 - ((H - K)) mag, (1) 
0' -H)> 2.5 - ((H - K)) mag. (2) 

As these pBzK galaxies have maximally-old, passively 
evolving systems, we derive a typical ((H — K)) color 
indicative of a stellar population with an instantaneous 
sta r-formation history forme d at Zf orm = 10 based on 
the lCharlot fc Bruzuall ([20071 hereafter CB07) models of 
((H - K))=0.7 mag. In Figure Q] we show the (B - z') 
and (V — H) colors with these criteria illustrated as a 
shaded polygon. To rule out any potential image arti- 
facts associated with the ERS field edges, we require the 
objects to be in the portion of the _ff-band mosaic which 
received the full complement of two orbits per pointing. 
We restrict the iJ-band magnitude to a relatively con- 
servative limit o f H < 25 mag to ensur e a complete and 
reliable sample (Wind horst et all 120101 ). With these re- 
quirements, we identify a final sample of 30 galaxies in 
43.1 arcmin 2 . Of these galaxies, 15 objects have (B — z') 
colors redder than the plot limits and are not represented 
in Figure [TJ We tabulate the 30 candidates in Table Q] 
and show color stamps of each object in Figure [2] 

We show the source counts of BzH galaxies in Fig- 
ure The counts plateau around H ~ 22.5 mag, which 
is roughly 4 mag brighter than the formal ERS complete- 
ness limit ([Windhorst et al.ll2010|) . suggesting that the 
faint-end of their luminosity function is relatively flat or 
declining. At the bright-end (H<22 mag), our survey 
and modified pBzK s election method p roduces source 
counts consistent with lLane et a l. (2007), based on the 
UKIRT Infrared Deep Sky Survey (UKIDSS) Ultra Deep 
Survey (UDS). While the UKIDSS/UDS survey is much 
wider (~ 0.6 deg 2 ), the ERS data pushes > 2 AB mag- 
deeper, into a regime not routinely possible from the 
ground. 

4. ANALYSIS 

4.1. Surface Brightness Models 

To determine the rest frame optical morphologies of 
our galaxies, we model the two dimensional light distri- 
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Fig. 1. — The BzH object selection. The small dots represent all 
objects detected and measured by SExtractor. Th e g rey area is the 
color selection region denned by equations [l] and l2lwith a typical 
color of ((H — K)) =0.7 mag. We identify 30 galaxies strictly meet- 
ing the BzH color criteria (shown as large filled circles), however 
15 galaxies have (B — z 1 ) colors too red to be shown here. There are 
several galaxies which formally meet our color criteria but were re- 
jected as being too faint or too close t o survey edges to be st udied. 
The r ed lines sho w the tracks for the | Colcman. Wu, & Wccdman 
l|l 980 . dashed) and lCharlot fc Bru zual (2001 solid) templates. We 
note, Galactic stars will span roughly similar (B — z 1 ) colors but 
will be >0.7 mag too blu e in (z 1 — H) to ha ve been misidentified 
as BzH galaxies (sec also Daddi et al. 2004). 



bution in the H-band using GalFit ([Peng et al.l 120021) . 
We determine an empirical point-spread function (PSF) 
from a median stack of 75 stars identified in the field, 
based on their full- width at half-maxim um (FWHM) and 
brigh tness determined by SExtractor ([Windhorst et all 
2010). We fit the standard Sersic profile: 



£(r) = E e e 



-b n [(r/r 



(3) 



where r e is the effective radiu^F 2 !. S e is the surface 
brightness at the effective radius, n is the Sersic in- 
dex, and b n is a constant found by numerically solv- 
ing r(2n) = 2"f(2n,b n ). Since GalFit is minimizing a 
goodness-of-fit, modeling the uncertainty in each pixel 
is critical for meaningful estimate of the uncertainty in 
each of the fit parameters. Therefore, we transform the 
weight maps produced by MultiDrizzle to uncertainty 
maps and explicitly include the shot noise of the objects. 
We perform all fits with GalFit in units of counts using 
the H-band mosaic effective exposure time of 5017.7 sec- 
onds. 

We excise a 81 pix x 81 pix stamp centered on each 
galaxy; this size was chosen as a compromise between a 
sufficient number of sky pixels for a robust sky estimation 
by GalFit, and the number of (generally unassociated) 
nearby galaxies. Neighboring galaxies can bias the pa- 
rameter estimation of the primary galaxy, therefore one 
must carefully mask out the unmodeled objects, or simul- 
taneously fit all the objects in the stamp. We opt for the 
latter, since it avoids the ambiguities in pixel masking, 
and permits the flux in a given pixel to be represented 
as the sum of several independent components. However 
without good initial conditions, the GalFit algorithm 

22 We refer to the effective radius as r e when in angular units 
and as R e when converted to physical units. 



TABLE 1 
BzH Sample and HST Photometry 
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408 


3 32 28.02 


-27 40 31.2 


27.52 ±0.49 


26.18 ±0.12 


25.18 ±0.08 


24.40 ± 0.05 


23.77 ±0.03 


22.87 ±0.01 


22.54 ±0.01 




606 


3 32 23.27 


-27 40 45.8 


> 28.11 


27.60 ±0.53 


26.12 ±0.23 


25.36 ±0.14 


25.16 ±0.17 


23.66 ±0.03 


23.10 ±0.02 




1696 


3 32 42.08 


-27 4141.2 


> 27.36 


> 27.29 


26.69 ±1.12 


26.16 ±0.77 


25.47 ±0.32 


24.98 ±0.14 


24.08 ±0.08 




2227 


3 32 42.34 


-27 42 04.0 


> 27.55 


26.85 ±0.63 


24.53 ±0.13 


24.50 ±0.14 


23.86 ±0.07 


22.91 ±0.02 


22.32 ±0.01 




2377 


3 32 25.03 


-27 42 09.6 


>28.10 


27.80 ±0.67 


27.68 ± 1.04 


27.43 ± 1.00 


26.44 ±0.42 


25.42 ±0.11 


24.52 ±0.06 X 


ray ID 226* 


2749 


3 32 14.92 


-27 42 21.9 


> 27.23 


26.40 ± 0.44 


24.92 ±0.19 


24.64 ±0.19 


24.67 ±0.19 


23.22 ±0.03 


22.69 ±0.03 




2750 


3 32 14.88 


-27 42 23.3 


> 27.26 


25.70 ±0.22 


24.56 ± 0.14 


23.92 ±0.09 


23.95 ± 0.09 


22.46 ±0.02 


21.97 ±0.01 X 


ray ID 145* 


2871 


3 32 31.09 


-27 42 26.6 


>28.23 


26.85 ±0.22 


25.78 ±0.14 


24.91 ±0.08 


24.17 ±0.05 


23.39 ±0.02 


23.07 ±0.02 




3000 


3 32 43.93 


-27 42 32.4 


26.23 ±0.47 


25.12 ±0.13 


23.61 ±0.06 


22.64 ±0.03 


21.99 ±0.02 


21.07 ±0.00 


20.59 ±0.00 




3152 


3 32 26.05 


-27 42 36.6 


> 27.96 


27.21 ±0.44 


24.93 ± 0.09 


23.93 ±0.05 


23.31 ±0.02 


22.15 ±0.01 


21.66 ±0.00 




3237 


3 32 35.91 


-27 42 40.9 


> 28.04 


26.77 ±0.30 


25.24 ±0.13 


24.60 ±0.08 


23.90 ±0.05 


22.67 ±0.01 


22.06 ±0.01 




3360 


3 32 30.58 


-27 42 43.4 


> 28.06 


27.27 ±0.44 


25.90 ± 0.22 


24.65 ±0.08 


24.11 ±0.05 


22.78 ±0.01 


22.22 ±0.01 




3376 


3 32 35.13 


-27 42 37.0 


> 28.04 


27.95 ±0.74 


26.80 ±0.45 


26.69 ± 0.48 


25.40 ±0.16 


24.91 ±0.07 


24.55 ±0.06 




3471 


3 32 27.94 


-27 42 45.7 


28.29 ± 1.34 


26.29 ±0.17 


24.75 ± 0.07 


24.02 ± 0.04 


23.48 ±0.02 


22.35 ±0.01 


21.91 ±0.00 




3488 


3 32 25.65 


-27 42 46.8 


> 27.97 


> 28.20 


> 27.63 


26.95 ±0.70 


26.01 ±0.31 


25.38 ±0.12 


24.93 ±0.10 




3551 


3 32 36.28 


-27 42 49.4 


27.78 ±1.42 


26.05 ±0.24 


24.32 ± 0.08 


23.38 ±0.04 


22.78 ±0.02 


21.63 ±0.00 


21.14 ±0.00 




3812 


3 32 36.66 


-27 42 58.5 


> 28.02 


26.05 ±0.15 


24.87 ± 0.09 


23.75 ± 0.04 


23.35 ±0.02 


22.23 ±0.01 


21.74 ±0.00 




4148 


3 32 44.97 


-27 43 09.1 


27.33 ±0.48 


25.00 ±0.05 


24.43 ± 0.05 


24.09 ± 0.04 


23.82 ±0.03 


22.63 ±0.01 


21.83 ±0.00 




4173 


3 32 41.24 


-27 43 09.7 


> 27.57 


26.31 ±0.30 


25.24 ±0.19 


24.48 ±0.11 


23.91 ±0.07 


23.29 ±0.03 


22.65 ±0.02 




4324 


3 32 31.32 


-27 43 16.1 


> 27.48 


26.23 ±0.29 


24.65 ±0.12 


23.85 ±0.07 


23.23 ± 0.04 


22.25 ±0.01 


21.77 ±0.01 




4327 


3 31 59.71 


-27 43 15.5 


> 28.40 


28.48 ± 1.29 


26.81 ± 0.47 


26.89 ±0.61 


25.90 ±0.36 


25.47 ±0.17 


24.37 ±0.08 




4534 


3 32 29.99 


-27 43 22.7 


> 28.07 


27.44 ±0.51 


25.42 ±0.13 


24.37 ±0.06 


23.84 ±0.04 


22.74 ±0.01 


22.26 ±0.01 




4648 


3 32 26.10 


-27 43 26.7 


28.39 ±1.28 


26.60 ±0.23 


25.47 ±0.13 


24.77 ±0.08 


24.29 ± 0.07 


23.49 ±0.02 


22.92 ±0.02 




4846 


3 32 25.93 


-27 43 31.1 


> 28.28 


> 28.37 


26.45 ± 0.31 


26.23 ±0.31 


25.28 ±0.17 


24.08 ± 0.04 


23.67 ±0.03 




4921 


3 32 02.44 


-27 43 35.8 


> 28.03 


> 28.14 


> 27.58 


26.98 ±0.75 


27.88 ±1.69 


26.13 ±0.23 


24.91 ±0.10 




5410 


3 32 32.11 


-27 43 55.3 


> 27.22 


26.33 ±0.38 


25.26 ± 0.24 


25.83 ±0.50 


25.53 ±0.43 


23.87 ±0.06 


23.45 ± 0.05 




5529 


3 32 01.77 


-27 44 01.1 


27.56 ±0.88 


28.02 ±1.16 


25.76 ± 0.24 


24.76 ±0.12 


24.39 ± 0.09 


23.46 ±0.03 


22.87 ±0.02 




5685 


3 32 12.78 


-27 44 07.7 


27.98 ±1.78 


27.70 ±1.10 


25.17 ±0.18 


24.95 ±0.18 


24.17 ±0.09 


23.44 ±0.03 


23.06 ± 0.03 




5735 


3 32 02.83 


-27 44 09.7 


> 28.05 


26.80 ±0.29 


26.88 ±0.52 


26.52 ±0.45 


26.10 ±0.33 


25.80 ±0.17 


24.66 ±0.08 




6845 


3 32 06.57 


-27 45 14.0 


> 27.57 


27.81 ±1.15 


25.93 ± 0.34 


25.06 ±0.18 


24.80 ±0.15 


23.44 ±0.03 


22.69 ±0.02 




7202 


3 32 06.40 


-27 45 54.7 


27.10 ±0.61 


25.76 ±0.15 


24.78 ±0.10 


23.67 ±0.05 


23.16 ±0.03 


22.37 ±0.01 


21.71 ±0.01 X- 


ray ID 82* 



t Coordinates refer to the J2000 epoch. 
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Fig. 2. — BzH color images. Each stamp is fs4''5 (~ 38 kpc at 2 = 1.6) on a side, has north up and east left, and has a pixel scale of 
0'.'090 pix" 1 . All images are shown with the same color and logarithmic intensity scale. We find the object just to the south of #2750 (and 
#2749) is just barely too blue to have been included in our sample, which is consistent with the suggestion that this object has recently 
(~150 Myr) undergone an intense burst of star formation of 500 — 2000 Mq yr _1 (van Dokkum & Brammcr 2010). 
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Fig. 3. — BzH galaxy counts. We show the differential surface 
den sity of ga laxies selected based on the usual BzK color crite- 
ria (Daddi ct al. 2004) and a typical color of {(H — K)) =0.7 mag, 
assuming an instantaneous burst of star formation and passive evo- 



lution from 2f 



■ 10. For comparison, we show the BzK object 



counts from La ne et al.l (120071) from the UKIDSS/UDS survey as 
open symbols, which were derived from a strict BzK selection. The 
upper axis is simply a transformation using the assumed {(H — K)) 
color and the right axis applies to our objects. To demonstrate that 
the apparent plateau in the BzH counts Hab ~ 22 mag is not due 
to simple survey completeness, we s how the total source c ounts 
from the ERS data as a dashed line (Windhorst ct al. 2010). The 
uncertainties presented here reflect the Poisson variation in the 
object counts, and do not include any contribution from cosmic 



variance. 



may not converge to a meaningful solution, since the 
number of degrees-of-freedom can be very large. There- 
fore, to estimate better initial conditions for each galaxy 
in the stamp and ultimately ensure convergence in the 
final simultaneous solution, we fit the two-dimensional 
light profiles of the primary BzH galaxy in question and 
any neighboring galaxies in a multi-stage process as fol- 
lows: 

1. We identify all neighboring galaxies and their 
associated pixels with SExtractor, using the 
same settings mentioned in § 12.21 Any sources 
with FWHM_IMAGE ff < 1 are eliminated from the 
SExtractor segmentation maps and object cata- 
logs, since they are likely not galaxies. 

2. The pixels of any galaxy whose isophotes are trun- 
cated (based on the SExtractor flags) are masked, 
by setting the uncertainty maps to 10 10 ADU. 
These galaxies are no longer considered in the 
GalFit process. 

3. We model the light distribution of each remain- 
ing galaxy (including the primary BzH galaxy) 
individually, while masking all the pixels associ- 
ated with every other galaxy. For any neighboring 
galaxy with semi-minor axis of B_IMAGE# < 1 pix 
from SExtractor, we switch from fitting a Sersic 
profile to a PSF model, in order to eliminate degen- 
erate degrees of freedom. Our results are robust to 
the choice of semi-minor axis limit, provided that 
we do not permit it to be larger than the size of the 
empirical PSF (discussed in more detail below). It 
is important to note that this step is only present to 



get reasonably accurate initial conditions for sub- 
sequent simultaneous object fitting. 

4. We refit a GalFit model, which contains a combi- 
nation of point sources and Sersic profiles, to the 
stamp as a whole using the results from the previ- 
ous step as the initial guesses. 

In a few cases, we manually masked diffraction spikes 
or stellar halos clearly associated with foreground stars, 
which were just outside the field-of-view of the stamp 
before proceeding through the above procedure. 

For the individual and simultaneous fits, we placed 
constraints on various GalFit parameters to the algo- 
rithm from diverging into an unphysical regime. We con- 
strained the centroid of any component to be within the 
±2<7 of the centroid determined by SExtractor, and the 
total magnitude to be 35 < H tot < MAG_AUT0 H - 2 mag. 
In general, these constraints are so weak that they gen- 
erally play no role in the fitting whatsoever, but they 
reduce the sensitivity to the pixel masking with the seg- 
mentation maps. We additionally constrain the Sersic 
index and effective radius to be 0.01 < n < 8 and 
0.01 < r e /A_IMAGE < 5, respectively. These constraints 
are considerably stronger, and we recognize that the 
model is likely incorrect when GalFit converges to a so- 
lution which is on these boundaries. These cases are rare, 
and generally a sign that additional astrophysical com- 
ponents are needed (e.g. bulge/disk separately, nuclear 
point sources, or merger signatures), that the frame was 
inappropriately sized, that the object was unresolved, 
and/or that there was some additional light component 
or image defect present in the image (e.g. diffraction 
spikes, stellar halos, and/or cosmic rays). 

As our primary interest here is on the sizes of these 
galaxies, it is imperative that we ensure the effective 
radii are robustly measured. While GalFit reliably 
determines the random uncertainty which follows from 
the maximum likelihood analysis, the total uncertainty 
should include a systematic term as well. To estimate 
the contribution from the systematic uncertainty, we 
construct a grid of simulated galaxies with brightnesses 
20 < H < 25 mag, effective radii 0.5 <r e < 5.5 pix, and a 
fixed Sersic index of n — A. These galaxies are convolved 
with the PSF and embedded in a blank region of the H- 
band mosaic. We then fit these simulated galaxies with 
GalFit using the above procedure, and find that the ef- 
fective radii are generally uncertain by ~ 10%, which is 
somewhat brightness dependent. We quadratically add 
this systematic uncertainty to the random uncertainty 
determined by GalFit. 

Many of the these red galaxies are very small and, even 
with the space based imaging, may still be unresolved. 
Therefore it is critical to properly identify which galax- 
ies are resolved and have reliable effective radii measure- 
ments. We begin by swapping the fully variable Sersic 
model for the primary galaxy with a pure PSF model, 
which can only vary in position and brightness. How- 
ever, as we have simultaneously fit every object in our 
postage stamps, we anticipate that GalFit may incor- 
rectly change the parameters of neighboring (unrelated) 
galaxies to compensate for the poor primary model, par- 
ticularly in the case of a well resolved primary galaxy. 
Therefore, when we use the PSF model for the primary 
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Fig. 4. — Comparison of PSF and Sersic models. We show the 
fractional difference between the x 2 goodness-of-fit measures from 
GalFit for the PSF and Sersic models. The BzH galaxies are 
shown as filled blue points, while the red asterisks represent 75 ERS 
point sources, presumed to be stars. For GalFit to optimally fit 
an unresolved source with a PSF, it will drive the effective radius 
and Sersic index to unphysical regimes, de facto fitting a PSF. 
Therefore, the goodness-of-fit of a Sersic model and a PSF should 
be roughly equal. However, if one were to fit a PSF to a resolved 
object, then there should be a noticeable increase in t he goo dncss- 
of-fit statist ic. W e adopt a critical value of -F cr it =0.01 ( Bond ct al. 
12003 ). see § l4,ll for more details. Of the 75 stars used in this test, 
only 4 have F > F CIit suggesting this is reliable way to classify 
unresolved objects. 

galaxy, we hold the parameters of the neighboring sources 
fixed at the values found previously by GalFit. We now 
compare the goodness-of-fit statistics for these Sersic and 
PSF models by considering the quantity: 



F 



X PSF x Sersic 



*Se 



(4) 



crsic 



and expect that sources with low values of F are equally 
well characterized by a PSF model as by the more com- 
plex Sersic profile. We calibrate this quantity by com- 
puting the .F-values of the known ERS stars, which were 
used to derive the empirical PSF used above. In Fig- 
ured! we show the F- values as a function of the GalFit- 
derived PSF magnitude, with the galaxies plotted as 
fillcd-blue points and the stars as red asterisks. The 
trend is as expected: nearly all Galactic stars can be 
found at — 0.5<-F<0, while the BzH galaxies are gen- 
erally at F > which depends on brightness. There- 
fore, we define objects which can be equally character- 
ized by a PSF as by a Sersic fit as having F < F crit , 
while objects with F > F cr it are more extended than the 
known ERS stars. We adopt F clit — 0.025 and note that 
only 4/75 s t ars ha ve F > F cr a , which is consistent with 
Bond et al.l (|2009fl . We give the GalFit results in Ta- 
bled 

4.2. Additional Photometry 

Our BzH criteria are designed to select passively 
evolving galaxies at z ~ 1.5, consequently the HST 
photometry covers only A ros t < 6000 A. We can ex- 
tend our wavelength coverage to A rcs t ~ 3 /jm with the 
K s -band imaging fr om the Very Large Telescope ( VLT, 
iRetzlaff et alj 12010) and the four IRAC channels from 
the Spitzer Space Telescope (Spitzer; PI: M. Dickinson). 



However the notably lower spatial resolution (0'.'7 — 2") 
of these images demands a different approach for mea- 
suring the flux, as source confusion can significantly bias 
the photometry. 

By inspection of the images, it is clear that our galax- 
ies are unresolved in the VLT and Spitzer data. There- 
fore, we obtain total magnitudes from a GalFit model 
with a similar approach discussed in § 14. li with a few 
simplifications. First, we assume that all the sources 
in a given postage stamp are unresolved, and are there- 
fore ideally modeled by a PSF. Second, we do not allow 
the centroids of all objects to vary more than ±1 pix (of 
the VLT or Spitzer images). Third, we only perform si- 
multaneous fitting with initial magnitudes given by the 
-ff-band measurements, as in many cases multiple HST 
sources are blended into a single VLT or Spitzer source. 

4.3. Stellar Populations and Photometric Redshifts 

We fit the 15-band (HST, VLT, and Spitzer) photom- 
etry with a library of stellar population synthesis mod- 
els to simultaneously determine the stellar mass, pop- 
ulation age, and redshift with our own software. Our 
model grid consists of a four dimensional parameter space 
spanned by redshift (z), stellar population agj^T (t), 
y-band extinction (Ay), and star formation timescale 
(r) for an exponentially declining star formation his- 
tory. We assume solar mctallicity, a Salpeter initial 
mass function, and adopt the CB07 population synthe- 
sis models. The allowed photometric redshifts ranged 
fro m < z < 7 with Az = - 01, th e ages adopted 
by iBolzonella. Miralles. & Pellol (|2000D . extinctions of 
< Ay < 2 mag with AAy =0.2 mag, and star formation 
timescales of r e [10~ 3 , 1,2,3, 5, 15, 30, 10 3 ] Gyr. We in- 
clude a systematic uncertainty on the observed fluxes of 
the form aj — axf, where we adopt a = 0.1, 0.1, 0.2, and 
0.2 for the ACS, WFC3, VLT, and Spitzer data, respec- 
tively to account for uncertainties in the zeropoints, tem- 
plates, and measurement approach (whether SExtractor 
or GalFit, for the HST or VLT/ Spitzer data). We com- 
pute the lcr uncertainties on the stellar population pa- 
rameters (e.g. mass, age, etc.) by use of a simple Monte 
Carlo calculation. For each band for a given galaxy, we 
draw a normal random variable with mean and standard 
deviation equal to the flux and flux uncertainty (with- 
out the the systematic term). By repeating for many 
iterations, we build up a distribution of each stellar pop- 
ulation parameter and take the mean and standard de- 
viation as the measured quantity and lcr uncertainties, 
respectively. 

We apply this approach to the complete sample of 
30 BzH galaxies. In Figure EJ we show the photometric 
redshift distribution and the comparison to a spectro- 
scopic redshift, when available. As expected, the me- 
dian redshift of the sample is z p hot = l-6±0.6, where the 
uncertainty reflects the standard deviation of the dis- 
tribution. To estimate our typical redshift uncertainty, 
we compare to the published spectroscopic redshifts in 
the lower panel of Figure Based on the root-mean 
squared scatter of the seven objects with known spectro- 
scopic redshifts, we estimate our uncertainties in (1 + z) 
are ~4.6%. We do recover the known galaxy cluster at 

23 We impose the usual self-consistency constraint that the a 
galaxy be younger than the age of the Universe. 
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TABLE 2 
BzH GalFit Results 



~W 



Notes 



408 


0.04 ± 0.00 


4.72 ± 0.45 


0.487 


0.160 






606 


0.09 ±0.01 


6.92 ±0.92 


0.404 


0.158 






1696 


0.41 ±0.02 


1.21 ±0.10 


0.356 


0.170 






2227 


0.35 ±0.01 


2.93 ±0.09 


0.364 


2.517 






2377 


0.16 ±0.01 


4.17 ±1.01 


0.322 


0.049 






2749 


0.51 ±0.02 


0.57 ±0.05 


1.849 


0.191 






2750 


1.03 ±0.19 


5.65 ±0.58 


1.891 


0.591 






2871 


0.03 ± 0.00 


8.00 ±2.76 


0.567 


0.023 


unresolved 


3000 


0.60 ±0.01 


8.00 ±0.14 


0.814 


9.440 


tidal tail 


3152 


0.14 ±0.00 


4.66 ±0.12 


0.515 


3.323 




3237 


0.22 ±0.01 


7.54 ±0.30 


0.432 


2.268 




3360 


0.11 ±0.00 


7.58 ±0.46 


0.466 


1.382 




3376 


0.15 ±0.01 


2.21 ±0.44 


0.349 


0.054 




3471 


0.08 ± 0.00 


4.40 ±0.16 


0.478 


1.525 




3551 


0.44 ±0.01 


2.74 ±0.11 


0.614 


0.851 


additional nuclear point soi 


3812 


0.16 ±0.00 


3.12 ±0.07 


0.486 


3.767 




4148 


0.04 ± 0.00 


5.73 ±0.32 


0.443 


0.690 




4173 


0.37 ±0.00 


0.98 ±0.02 


0.359 


2.373 




4324 


0.25 ±0.00 


3.27 ±0.08 


0.404 


4.536 




4327 


0.03 ±0.02 


2.96 ±3.62 


0.411 


0.000 


unresolved 


4534 


0.14 ±0.00 


2.21 ±0.07 


0.388 


2.718 




4648 


0.17 ±0.00 


2.08 ±0.10 


0.355 


0.913 




4846 


0.16 ±0.00 


1.30 ±0.13 


0.506 


0.232 




4921 


0.19 ±0.08 


8.00 ±4.50 


0.337 


0.012 


unresolved 


5410 


0.00 ±1.00 


2.64 ±1.00 


0.364 


-0.009 


likely three distinct clumps 


5529 


0.29 ±0.00 


2.63 ±0.12 


0.369 


1.107 




5685 


0.48 ± 0.00 


0.36 ±0.02 


0.386 


1.399 




5735 


0.08 ±0.01 


8.00 ±4.15 


0.349 


0.012 


unresolved 


6845 


0.42 ±0.01 


1.33 ±0.04 


0.385 


1.648 




7202 


0.30 ± 0.00 


2.29 ± 0.04 


0.434 


5.825 







^The Sersic index in equation ((3j - 

J The fractional difference between the goodness-of-fit for the PSF and Sersic models (see S I4.1l for more details). 
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Fig. 5. — Photometric redshift distribution. In the upper panel, 
we show the photom etric redshift distribution derived by the pro- 
cedure described in § 14.31 As expected, our sample of BzH galaxies 
is generally located at z p hot = 1-6. In the lower panel, we show the 
fractional difference between our redshift estimates and spectro- 
scopic redshifts, where available. Based on the RMS on the frac- 
tional differences, we estimate that (1 + z) is accurate to ~4.6%. 



z ~ 1.6 ()Kurk et al.ll2009l ). We present our resulting pho- 
tometric redshifts and stellar population parameters in 
Table [3] 



5. BZH GALAXY SIZE EVOLUTION 

Based on the photometric redshift estimates and the 
spectroscopic data (where available), the BzH selection 
reliably identifies galaxies in the interval (z) ~ 1.6 ± 0.6. 
However to study the evolution of their sizes with red- 
shift, we must compare to other similarly selected sam- 
ples. The high redshift (z > 1.5) samples are generally 
derived from similar color criteria presented here, and 
have effective radii measured in the _ff-band, for a rest 
frame wavelength of A res t ~ 6500 A. To ensure fair com- 
parisons with lower redshift samples, we require simi- 
lar rest frame sizes. For the low redshift data, we use 
the sample of 8666 early-type galaxies at z ~ 0.2 with 
effective radii measured in the i'-band (|Bernardi et al.l 
2003) selected from the Slo an Digital Sky Su rvey, Early 
Data Release fSDSS-EDR: lYork et al.l[20l0l) . The stel- 
lar masses for the 7th Data Re l eas^ l galaxies were de- 
termined following Sali m et al.l ([2007). By cross match- 
ing these samples, we obtain 8595 galaxies for our low- 
redshift comparison. We select three mass ranges which 
are volume-limited based on the upper and lower flux 
limits imposed by de tector saturation in SDSS and the 
iBernardi et al.l (|2003f) brightness criterion. In Figure [5] 
we show these volume-limited selections (black boxes) 
and the flux limits (dashed lines) . If these limits are not 

24 Obtained from http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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TABLE 3 

Photometric Redshifts and Stellar Population 

Parameters 



ID 



^phot 



Age 1 " 
(Gyr) 



All 
(10 11 M ) 



A u 



408 

606 

1696 

2227 

2377 

2749 

2750 

2871 

3000 

3152 

3237 

3360 

3376 

3471 

3551 

3812 

4148 

4173 

4324 

4327 

4534 

4648 

4846 

4921 

5410 

5529 

5685 

5735 

6845 

7202 



1.59 



+U.U5 
-0.04 
i o 7 +0.24 
1 - S '-0.20 
4 12+ 056 
i 97 +0.11 
1 '^'-0.13 
2 S6+ - 13 
, ocr+0.27 
1 -°°-0.09 

l- 8 '-0.09 
1 54+ - 08 

1 T5+ - 42 

i no + 0.21 

■"-^O—n.ip 

1 66+ - 08 
1 - DD -0.06 

1 63+ - 44 
1 -°°-0.09 

1 49+ - 09 
1 -* a -0.11 

1 67+ - 06 

1 60+ - 06 
i - DU -0.10 
i 77 +0.07 
'•■' '-0.07 
9 OO+0.07 
z -°°-0.05 

0.98i°;i 7 

i E5+O.IO 
J - oo -0.09 
n 1K+1.10 
o - lo -0.77 

1 51+ - 17 
J - O1 -0.26 
n41 +0.61 

u - 4i -o.io 

o 7 +1.38 

°-°' -1.69 
9 i +0.17 
z ' lz -0.37 
, n o+0.24 
l-l«_o.35 
,+0.13 
-0.14 

2 64+ - 31 
Z - D ^-0.15 

1 67+ ' 64 
1 -°'-0.37 
i i 7 +0.11 



1.41" 



367 

615 



610 



611 



6111 



329 



c+U.U 
°-0.0 
7+0.3 
'-0.2 
+3.1 

0.2 
+ 1.0 

1.0 
+3.0 

0.4 
+0.0 

0.0 
c+0.2 
°-0.4 
7+0.0 
'-0.4 
+0.0 

5.0 
+0.0 

3.5 
+0.0 

0.2 
+ 1.0 

0.5 
+0.5 

0.1 
+0.0 

0.3 
+0.4 

0.3 
+0.0 

0.0 
+0.0 

0.0 
+6.0 

2.1 
+0.4 

0.3 

2.1 

0.1 
+3.1 

0.4 

5 +i.o 



7.1 

+4.0 
0.3 

2.8 

0.5 

+3.1 

0.3 

+3.0 

3.5 

1.2 

0.4 

>+2.0 

'-0.2 

5+1.9 

'-1.9 

;+2.9 

'-0.3 



0.4 + 0.0 
0.4 + 0.0 
11.3 + 1.1 
1.3 + 0.1 
0.1 + 0.0 
0.8 + 0.0 
0.8 + 0.0 
0.2 + 0.0 
5.1 + 0.2 
2.7 + 0.1 

1.1 + 0.1 

1.2 + 0.1 
0.1 + 0.0 
0.7 + 0.0 

2.2 + 0.1 
1.0 + 0.0 

1.3 + 0.1 
0.3 + 0.0 
0.9 + 0.0 
1.9 + 0.2 
0.9 + 0.0 
0.0 + 0.0 
0.3 + 0.0 
2.3 + 0.3 
0.2 + 0.0 
0.4 + 0.0 
0.2 + 0.0 
0.1 + 0.0 
0.9 + 0.1 
1.7 + 0.1 



1.5 
3.2 
2.1 
5.0 
0.6 
6.0 
8,1 
1.9 
4.1 
0.3 
1.0 
1,-1 
4.3 
2.5 
0.3 
4.9 
2.8 
1.8 
2.5 
2.3 
2.8 
2.8 
7.1 
2.0 
7.9 
3.5 
9.3 
3.0 
5.2 
2.0 



t Population age assuming an exponetial star formation history. 
* Stellar mass. 

strictly imposed, then an artificial redshift dependence 
on the effective radii will be introduced as the radii are 
tightly correlated with the stellar masses, which roughly 
scale with luminosity. For example, without these vol- 
ume limits, the effective radii of the SDSS galaxies will 
seem to increase with redshift since the survey is not sen- 
sitive to the lower mass (smaller) galaxies at the higher 
redshifts. Eliminating this potential Malmquist-like bias 
in the low redshift sample ensures a fair comparison be- 
tween the high and low redshift data. 

To investigate the passively evolving galaxy size evo- 
lution at a fixed stellar mass, we show in Figure [7] the 
effective radii as a function of redshift for this work (large 
circles ), the SDSS s ample (small dot s) . ILonghetti et all 
(120071 diamond s ). iDamianov et al.l (120091 t riang les) , 
iRutkowski et al.l (120101 crosses). iDaddi et all ([20051 as- 
terisks). and lCimatti et all (|2008l . squares). We overplot 
the two canonical models of (1 + z)~ a (solid line) and 
H{z)^ 13 (dot-dashed line), where H{z) is the Hubble pa- 
rameter, and present the best fit parameters in Table |3J 
While the data do not significantly favor either model, 
they serve to highlight an important trend: the amount 
by which galaxies are smaller in the past depends on the 
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Fig. 6. — Selection of SDSS galaxy sample. We show the stel- 
lar mass as a func tion of sp ectroscopic redshift for the early-type 
galaxy sample of Bcrnardi ct al. (2003) drawn from the SDSS- 
EDR. The stellar mass estimates wer e derived f rom stellar pop- 
ulation fits to the SDSS photometry (Salirn et al. 2007), and are 
similar in nature to those described in ^ 14.31 The dashed lines indi- 
cate the bright (i' < 14 mag) and faint ( i' > 17 mag) completen ess 
limits set by detector saturation and the Bcrnardi ct al. (2003) se- 
lection, respectively. The boxes show our volume-limited selections 
for proper low-redshift comparisons. If such limits are not imposed, 
then an artificial trend in effective radius with redshift will arise 
from a Malmquist-type bias. 



TABLE 4 
Mass-Dependent Size Evolution Models 



log A/ low 
(logM Q ) 



log M high 
(logM ) 



W 



Re(z = 0) 

(kpc) 



10.3 
10.7 
11.1 



10.7 
11.1 
11.5 



0.38 
0.87 
1.42 



0.49 
1.07 
1.86 



2.76 
4.36 
7.60 



tFor the (1 + z)~ a model. 
*For the H(z)- 13 model. 

stellar mass. iNewman et al.l ()2010[ ) identified a similar 
result for 12 galaxies with dynamical masses determined 
from measured velocity dispersions. 

Our stellar mass-dependent size evol ution model is 
qualit atively similar to that proposed by Hopkins et al. 
(2009), where the power-law index varies with mass as 
a{M*) sa 0.23 log (M»/10 9 M Q ). However, their model 
underpredicts our measured power-law indices in our 
three mass regimes, but does give the same qualitative 
steepening of a with mass. Following their approach, we 
derive tentative relationships for the power-law indices: 



a(M*)«-1.8+1.4x log 



/8(M.) w -2.3 + 1.8 x log 



AL 



10 9 M Q 

( M. 

Vio 9 m q 



(5) 



(6) 



and caution that with only three independent mass bins, 
these fits should be considered preliminary at best. How- 
ever, these resul ts give indices c onsist ent with the re- 
ported value of iBuitrago et al.l ()2006[ ). who find that 
a(M*> 10" M )wl.5. 



6. DISCUSSION 
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Fig. 7. — BzH galaxy size evolution. We show the effective 
radius versus redshift for the various stellar mass slices described 
in Figure [6] and several early-type galaxy samples. We do not 
show the uncertainties for clarity. As solid and dashed lines we 
show the best fit model of R e oc (1 + z)~ a and R e oc H(z)~P , 
respectively. While the data cannot rule out either model, we show 
both to help illustrate the emerging trend: the increase in effective 
radii from z ~ 2.5 to present is strongest in high mass galaxies 
(M* > 10 11 Mq). Based on these fits, we give a tentative estimate 
for a(M,) and /3(M«) in equations J5} and J5J, respectively. 

It has become relatively well established that early- 
type galaxies are indeed smaller at high redshift than 
their local co unterparts for a given stellar or dynam i- 
cal mass (e.g. iTruiillo et aHl2009l iNewman et al.N2010D . 
However the causes for this result are far less clear, or 
agreed upon. The likely mechanisms can be broadly 
characterized as astrophysical effects (early time major 
mergers, late time minor mergers, adiabatic expansion, 
or stellar mass-to-light gradients) , or observational biases 
(underestimating the effective radii, overestimating the 
stellar masses, or incorrectly assuming that the high and 
low redshift populations are directly comparab le). Using 
a semi-analytical model, Hopkin s et al.l (|2010a^ find that 
the factor of ~ 5 increase in the effective radii over the 
last ~ 10 Gyr for galaxies with M* > 10 11 M Q can be 
explained by a combination these effects, with the late- 
time minor merging playing the largest role. However, 
we find that lower mass systems exhibit notably weaker 
redshift evolution, suggesting a different mixture of the 
main processes may be at work. For example, the adi- 
abatic expansion mode may become more critical given 
the shallower gravitational potentials in the low mass 
systems. 

As noted above, the stellar mass- depen dent size evolu- 
tion se en here is similar to the findings of INewman et all 
(2010) for dynamical masses estimated from velocity 
dispersions, and has an interesting conse quence for 
the R e - M* relation. iShen et all (|2003t ) find that 
SDSS elliptical galaxies follow the scaling relation R e = 
i?e,ii(M*/10 n M q )t, where R e>n =4.16 kpc and 7 = 
0.56. However the mass-dependent size evolution ob- 



served here and by INewman et al.l ([20101 , imply either 
a fundamental change in the R e — M scaling relationship 
at high redshift (such as a flattening of the effective ra- 
dius for decreasing stellar mass) or a redshift dependent 
value of 7, in addition to the usual lower value of R e ,n- 
However, at this stage the high redshift data cannot dis- 
tinguish these two scenarios or shed light on the cause 
for this change. 

The red galaxy formation paradigm wherein merg- 
ers of gas-rich discs trigger intense starbursts which are 
later quenched by active galactic nuclei, finally giving 
way to dead sphe roidal systems spent of their gas (e.g. 
iFaber et al.l 120071 and references therein) suggests the 
passive galaxies may have come from a population of 
more active galaxies in their recent past. The Ly man- 
break galaxies (LBGs) observed at z > 3 are a possi- 
ble progenitor system. While their stellar masses are 
far less well constrained due to the lack of rest frame 
IR data, the majority of LBGs have M* ~ 10 10 M w be- 
tween 3 < z < 6 (e.g. [Papovich. Dickinson, fc Ferguson! 
[20011: 1 Yan etaLl [20061) . Therefore, the typical LBG at 
3 <z< 6 belongs in the lower panel of Figure [71 an d have 
effective radii of R e ~ 1 kpc (e.g. iFerguson et al.1 2004; 
Bouwens etlD 120061 : lHathi. Malhotra. fc Rhoadsl 120081 : 
Oesch et aLll2010f> . afte r transforming to a n equivalent 
rest frame wavelength ([Barden et al.l 120051) . The LBGs 
are then consistent with our H{z)~@ model, possibly sug- 
gesting that the physical mechanism driving the passive 
galaxy evolution may be at work for the LBGs as well. 
Given their increased star formation rates, larger gas con- 
tent, and lower tota l mass, the mass -loss modes, whether 
driven by AGN (e.g.lFan et al.ll2008l) or stellar winds (e.g. 
iDamianov et al.ll2009D . are likely more important. How- 
ever, we recognize that these samples (LBGs and passive 
galaxies) are quite different and that the next generation 
of space-based infrared instruments (Near- Infrared Cam- 
era and Mid-Infrared Instrument on James Webb Space 
Telescope) will provide a much clearer picture for the 
high redshift (z>3) size comparisons. 

7. SUMMARY 

We identified 30 passively evolving galaxies to H < 
25 mag from a set of color criteria similar to the -pBzK 
selection. We measure rest frame optical (A res t ~6500 A) 
effective radii as two-dimensional fits to the iJ-band im- 
age. By comparing with several other comparable sam- 
ples at various redshifts, we find that the size evolution 
depends on the stellar mass. We give tentative scalings 
between the power-law index of the R e — z relation and 
stellar mass. Future surveys, such as the coming Multi- 
Cycle Treasury programs with HST, will have the unique 
opportunity to improve upon our scalings and extend to 
both higher redshifts and lower stellar masses. 
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